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Abstract 

The well known (3+1) decomposition of Thorne and Macdonald is invoked 
■ to write down the Einstein-Maxwell equations generalised to (d+1) dimensions 

O 

bjo: 



and also to formulate the plasma equations in a flat FRW like spacetime in 
higher dimensions (HD). Assuming an equation of state for the background 
metric we find solutions as also dispersion relations in different regimes of the 
, universe in a unified manner both for magnetised(un) cold plasma. We find 

that for a free photon in expanding background we get maximum redshift in 
4D spacetime, while for a particular dimension it is so in pre recombination 
era. Further wave propagation in magnetised plasma is possible for a restricted 
q-i frequency range only, depending on the number of dimensions. Relevant to 

5^ ' point out that unlike the special relativistic result this allowed range evolves 

with time. Interestingly the dielectric constant of the plasma media remains 
constant, not sharing the expansion of the background, which generalises a 
' similar 4D result of Holcomb-Tajima in radiation background to the case of 

• higher dimensions with cosmic matter obeying an equation of state . Further, 

Tj- , analogous to the flat space static case we observe the phenomenon of Faraday 

rotation in higher dimensional case also. 
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1. Introduction 

While great strides have been made by general relativists to address the issues com- 
ing out of the recent observations in the field of astrophysics and cosmology and 
despite the fact that more than 90 percent of the cosmic stuff in stellar interior and 
intergalactic spaces is made up of matter in plasma state the much sought after 
union between the plasma dynamics and general relativity still remains elusive. Al- 
though we occasionally come across stray works like 'plasma suppression of large 
scale structure formation in the universe' [1] as also even the simulation techniques 
[2] where it is argued that the geometry of the spacetime can be purported to be 
structured by observing sound waves in primordial plasmas the fact remains that 
as the field equations both in general relativity and plasma dynamics are highly 
nonlinear it is very difficult to obtain closed form solutions in physics. So either 
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numerical relativity or a linearised approximation of the plasma equations is pre- 
ferred, before going in for more complicated non linear phenomena. If we briefly 
trace the thermal history of the universe theorists believe that from approximately 
t = 10~ 3 s to t — 1 s at temperatures T > 10 10 K there was an electron positron 
plasma at ultra relativistic temperatures. With cooling the plasma comprises mainly 
electron and hydrogen ions, may be with small amount of helium and other light 
elements in thermal equilibrium with photons. This period has come to be known 
as radiation dominated era because it is believed that the energy density of photons 
much exceed that of matter [3] . As the temperature decreased further with expan- 
sion at around t ~ 10 13 s recombination of ionized hydrogen atoms took place with 
consequent decoupling of matter and radiation and at a certain stage the universe 
becomes matter dominated. Consequently the temporal behaviour of the scale fac- 
tor correspondingly changes with the evolution of matter at different time scales. 
Nevertheless, the studies of the effects of this expansion on the electromagnetic in- 
teractions of the matter, including the longitudinal and transverse modes have not, 
so far got the legitimate attention it deserves. 

Another area of current interest is the role of an external magnetic field in many 
astrophysical systems and cosmology and possible sources for the field in different 
scales [1] although the origin of the field continues to evade plausible physical ex- 
planations [5] so far. Granted that temperature (particularly at early universe) is 
a known enemy of magnetism there are no compelling reasons why magnetic fields 
should not have been present in the early Universe either. Indeed, the presence of 
large scales magnetic fields in our observed Universe is a well established experimen- 
tal fact. Since their first evidence in diffuse astrophysical plasmas beyond the solar 
corona [4, 5] magnetic fields have been detected in our galaxy and in our local group 
through Zeeman splitting and through Faraday rotation measurements of linearly 
polarized radio waves. The Milky Way possesses a magnetic field whose strength 
is of the order of the microgauss corresponding to an energy density roughly com- 
parable with the energy density today stored in the Cosmic Microwave Background 
Radiation (CMBR) energy spectrum peaked around a frequency of 30 GHz. Faraday 
rotation measurements of radio waves from extra-galactic sources also suggest that 
various spiral galaxies are endowed with magnetic fields whose intensities are of the 
same order as that of the Milky Way. The existence of magnetic fields at even larger 
scales (intergalactic scale, present horizon scale, etc.) cannot be excluded, but it 
is still quite debatable since, in principle, dispersion measurements (which estimate 
the electron density along the line of sight) cannot be applied in the intergalactic 
medium because of the absence of pulsar signals. Given the fact that a seed field 
does exist its amplification mechanism through the so called dynamo effect is rel- 
atively well understood [6] for an expanding cosmological model with the help of 
magnetohydrodynamical (MHD) equations [7]. On the other hand when discussing 
a primordial magnetic field one should also bear in mind that a large value would 
create significant anisotropy of the background geometry [8] with consequent im- 
pact on CMBR findings. As the large scale isotropy of the CMBR fairly excludes 
that possibility there should be efficient mechanism to rapidly damp that field. 

On the other hand there has been, of late, a resurgence of interests in physics in 
higher dimensional spacetime [9] in its attempts to unify all the forces in nature, to 
give a physical explanation of the current accelerating era of the universe without 



bringing in any hypothetical quintessential type of scalar field [10] by hand, in the 
newly fashionable area of brane cosmology [1 1] where the gravity is supposed to act 
in the bulk while other forces in the physical 3D space. It has also received serious 
attention in the recently proposed induced matter theory pioneered by Wesson and 
others [12]. Most importantly both higher dimensional spacetime and cosmological 
plasma have one thing in common - both are very relevant in the context of early 
universe. While early universe may be loosely viewed as the history of evolution 
of matter in plasma state it can also be shown that starting from a higher dimen- 
sional phase the Einstein's generalized field equations dictate results such that as 
time evolves the 3D space expands while the extra dimensions shrink till plackian 
length when some stabilizing mechanism ( for example, quantum gravity, casimir 
effect, a repulsive potential) halts the shrinkage down to a very small length, say 
planckian size as to be invisible with the low energy physics at the moment. So the 
world around us appears manifestly three dimensional. But it should be emphasized 
that the time scales for the spontaneous self compactification (SSC) and the onset 
of nucleosynthesis clearly differ with SSC occurring much earlier. So many of our 
findings in MHD lose much of their relevance when working in higher dimensional 
spacetime. While literature abounds with works on the effects of the expanding 
background on the matter distribution and vice versa as also on a large number of 
other physical processes scant attention has been paid so far to address the issues 
resulting from the expanding universe on the propagation of say, electromagnetic 
wave as also its interactions in a plasma media. To authors' knowledge Holcomb 
and Tajima (HT) [13, 14], Banerjee et al [15] and later Dettmann et al [16] made 
important contributions in this regard. HT investigated the electromagnetic wave 
propagation in a radiation dominated and later in a matter dominated background 
both with or without any plasma material and also generalised it to the magnetized 
case both warm and cold. For the ultra relativistic case it is observed that all the 
modes redshift at the same rate i.e., photons are, so to say, self similar. Though not 
explicitly pointed out in their works we think this, however, may be a direct conse- 
quence of the conformal flatness of the FRW metric chosen by them. Later Banerjee 
et al discussed this in a little more general way. Dettmann et al got the similar results 
starting from a kinetic theory approach. HT's findings for the matter dominated 
model, however, differ sharply from the first paper in the sense that while in the 
unmagnetised plasma case the photons redshift identically but here for the Alfven 
waves the frequency redshifts in a bizarre fashion unlike the case of free photons, 
depending on the magnitudes of the plasma density and also strength of the external 
magnetic field. On the other hand Dettmann shows that even for the unmagnetised 
plasma the plasma oscillations and the photons do not share the identical temporal 
dependence if they are decoupled. They, however, treated the whole situation from 
kinetic energy considerations. In the present work we have investigated the plasma 
dynamics in an expanding higher dimensional background in a very general way. 
Taking a (d+1) dimensional flat FRW type of metric as background, which one of 
us [17] derived earlier in a different context we first take the case of propagation 
of an electromagnetic wave in vacuum. To make things very general we assume an 
equation of state for the background as p = 7p. Taking 7 = ~, (as we are dealing 
with a (d+1) dimensional case) for the radiation and 7 = 0, for the dust case in the 
resulting solution we observe that solutions are very similar to the special relativistic 



form except that the frequencies red shift depending on 7 and the field amplitude is 
no longer a constant decaying with the expansion rate, reminiscent of the acoustic 
case where the damping occurs due to some form of dissipating mechanism. Here 
the expansion of the background, in a sense, takes the role of dissipation, causing 
this type of damping. It is observed that with number of dimensions the red shift 
decreases, being maximum in the usual 4D case. Moreover, red shift is less in the 
post recombination era compared to the early universe as expected. After briefly 
carrying out the so called (3 + 1) decomposition of the Maxwell's equations gener- 
alised to (d + 1) dimensional spacetime in section 2 we investigate the propagation 
of an electromagnetic wave in vacuum in section 3. In section 4 the electrostatic 
oscillations are very briefly discussed for a cold plasma and the dispersion relations 
obtained both for transverse and longitudinal modes. In section 5 we discuss in some 
detail the propagation of an electromagnetic wave in cold plasma in the presence 
of an ambient and homogeneous magnetic field both parallel and perpendicular to 
the wave propagator k. The presence of the magnetic field introduces newer and 
interesting modes of oscillation, creating Left and Right circularly polarized waves, 
resulting in the well known classical phenomenon of Faraday rotation. We also find 
that the wave propagation in the plasma is possible for certain range of frequencies 
only and this range critically depends on the number of dimensions. The paper ends 
with a short discussion in section 6. 



2. Field Equations and its Formalism 

We extend here the (3+1) decomposition of GR as formulated by Arnowitt, Deser 
and Misner (ADM) [18] to a higher dimensional space time of (d+1) dimensions. 
The ADM formalism was developed mainly to address the issue associated with 
numerical relativity as also quantization of gravity fields, specially when the space 
time has considerable symmetry. This work is connected with plasma physics in 
curved space time. To make use of the intuition from the known results of MHD 
in flat spacetime it is preferable to spilt the ordinary electromagnetic field tensor 
F^ u into electric and magnetic fields E and B in terms of which the field equations 
are more familiar. As the split formalism has been extensively discussed and used 
in the literature [13, 14, 19] we shall very briefly touch upon its salient features as 
extended to higher dimensions. 

We define a set of observers ( Fiducial Observers or FIDO ) at rest in the space 
spanned by the hypersurfaces of constant universal time, having a d-velocity vector 
field, n orthogonal to spatial slices. It is well known that for a rotation-free space 
time ( as we are dealing here ) 

= CTij + ^ - (l >".i (!) 

(i, j = 1,2,3, ... d) where cr^ is the shear of the Eulerian world lines given by 
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Here 7^ is d- metric spatial tensors 

a* = n%v? (3) 



and 

9 = nj(= -K) (4) 



are the usual acceleration and expansion scalars and K is the trace of the extrinsic 
curvature, K\. 

For our metric we take the ( d+1) dimensional generalized FRW space time as 
ds 2 = dt 2 - A 2 (dx 2 + dy 2 + dz 2 + difj 2 n ) (5) 

(n = 5, 6, 7, . . . , d ) 

where A = A(t) is the scale function. For this particular metric two relevant 
quantities a — % ( the lapse function - the rate of change of fiducial proper time 
to that of universal time ) and also the shift vector (3 ( signifying how much spatial 
co-ordinates are shifted as one moves from one hypersurface to the other ) naturally 
reduce to a = 1 and (3 = 0. 

In an earlier work [17] one of us extensively discussed the ( d+1) dimensional 
isotropic and homogeneous space time and assuming an equation of state, P = 7p 
found the scale factor as ( P = pressure, p = energy density ) 

A ~ t^i = t n , n = — 2 (6) 

d(l + 7) 

With the extrinsic curvature scalar defined as 

K = -9 = -d- A = -d- (7) 
At w 

we finally write down the Maxwell's equations [20, 21, 22] (see Mcdonald et al 
for ( 3+1 ) split for more details ) generalized to ( d+1) dimensions as 



V.£ = 4vrp e (8) 

V.B = (9) 
BE 

— = KE + cA^V x B -4?rJ (10) 
at 

Ft R 

— = KB — cA~ l V x E (11) 
dt v ; 

= Kp e — V.J (charge conservation) (12) 



and finally the particle equation of motion in (d + 1) dimensional as 

A d - X q (E + A^x B^j (13) 



DA d ~ l p 



Dt 
or 



Dp d — 1 / v 



Dt d 

where 



Kp + q^E + A^xB^j (14) 



= ~ (dt + v.V) (15) 
Dt a 



is the convective derivative and the d- momentum 



p = m e Tv (16) 

( m e is the rest mass, Y is the boost factor, and v, the d- velocity ). We thus see 
that for the simple metric given by equation(5)the decomposed ( d+1)- dimensional 
Maxwell equations closely mimic the flat space counterparts with some additional 
inputs from curved geometry( e.g., A and K terms). 

Here V- and Vx are the ordinary Minkowskian divergence and curl in Cartesian 
co-ordinates. In what follows we shall consider, for simplicity, the small amplitude 
linear theory such that the convective derivative simply reduces to ordinary deriva- 
tive, | . 



3. Electromagnetic Waves in Vacuum 

With the set of equations split to ( d + 1) formalism we are now in a position 
to attempt applications in varied plasma phenomena. To start with we discuss 
briefly the propagation of an electromagnetic wave in free space in the expanding 
back ground. Using equations (9-11) we get via V • E — (for vacuum) the wave 
equation as 

A 2 E + (2d + 1) AAE + (d 2 A 2 + dAA\ E = c 2 V 2 E (17) 
Assuming separation of variables in electric field 



£(x, t) = E t E r 



we finally get 



E, 



(2d+l)jE t 



A 2 



+ 



E t = 



(19) 



which, through equation (6) finally reduces to 



t 2 E t + (2d + l)ntE t + [dn (dn + n-l) + k 2 c 2 t 2{1 ~ n) ] E t = (20) 

(here ki is a separation constant) corresponding to some initial fiducial time tj. 
We shall subsequently see that k^ is also identified with the d- dimensional wave 
vector. 

On the other hand, as we are dealing with a homogeneous world the spatial equa- 
tion remains unchanged yielding a solution e l<ykl ' r -* as in special theory of relativity. 
A little algebra shows that the time equation is reducible to a Bessel equation of or- 
der | as in the 4D case. Thus dimensionality or the equation of state has apparently 
no role in determining the order of the equation. Plugging everything together we 
get 



l-(2d+l)n (r,\ 

E = E et 5 nr 

2 



1.1— n 



1 — n 



-kiC 



iki.r 



(21) 



Hi we get 



(2) 1 

where H\ is a Hankel function of order ^. Replacing the asymptotic form of 
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_ „ / 2kicd(l + 7) 2 ifc t cd(i+ 7 ) , d( d 1 /|+^ 2 

V w{d(l + 7) - 2} 1 J 

It represents a d-dimensional 'damped' harmonic wave as one encounters in me- 
chanical vibration. While in mechanical motion the damping occurs due to friction 
here the expansion of the universe seemingly causes some sort of damping. For 
pre recombination era in (d + 1) dimensional space time 7 = \ and the damping 
factor is A~ d . Hence the damping of the wave amplitude apparently decreases with 
number of dimensions. This finding merits some explanation. We have remarked 
earlier that the amplitude decay is somewhat geometrical in nature caused by the 
expansion and curvature. In that with dimension the expansion rate de- 

creases one expects that damping should be larger in 4D but a little inspection of 
the last relation shows that when we plug in the expression of 'A' (equation (6)) the 

2d 

last relation further reduces to t d + x . So the damping actually increases in higher 
dimensional spacetime. On the other hand, for a fixed d the damping factor is t~ 2 
for post recombination era (7 = 0) (see figure 1). Alternatively for the case of a very 
large number of dimensions the damping asymptotically reaches t~ 2 , a form set for 
post recombination era. It has not escaped our notice that the scaling of E or B 
for 7 = is independent of the number of dimensions unlike the radiation case. So 
the amplitude factor gets increasingly damped as the universe ages. The fact that 
damping of the E or B increases with the number of dimensions in the early universe 
has a number of interesting theoretical implications. Firstly we mention in the intro- 
duction that in order that the universe evolves isotropically according to the FRW 
model there should be efficient mechanisms to damp the primordial magnetic field 
as early as possible. In that respect the HD spacetime has some inherent advantage 
over the standard 4D in the sense that the damping is faster in HD. Secondly the 
magnetic field at small scales may influence the bigbang nucleosynthesis and change 
the primordial abundances of light elements by significantly changing the expansion 
rate of the universe at the corresponding time. The success of the standard BBN 
scenario can provide an interesting set of bounds on the intensity of the magnetic 
field at that epoch [5], indirectly constraining the number of dimensions of the 
spacetime. At this stage it may not be out of place to call attention to the fact that 
most of the above findings are of theoretical nature only and it is not feasible to 
relate them to current astrophysical data. Because multidimensional cosmological 
models lose much of their relevance well before the onset of bigbang nucleosynthesis 
and current observational findings can be explained for all practical purposes if the 
cosmological evolution be modelled along the standard four dimensional spacetime. 

If as usual we set kiC = uj.i ( the angular frequency of the wave at some initial 
time t = ti ) then the above equation may be rewritten as 



E = E n U I r y (1 ± 7) , r^e-^^ rTOJ e^ 
n{d(l + 7) - 2} 



^^0^f^-^'^ (23) 
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Figure 1: | Ed |~ t graph 



where 

u d = u i t~' s ^+^ (24) 
gives a measure of the red shift of the photon due to background expansion. For 

1 2 

radiation dominated era 7 = -j , u d = LUjt d +! , so the rate at which the frequency 
decreases is maximum in 4D universe. Moreover damping is greater in radiation era 
(see figure 2). 

Returning again to the equation (22) we see that the horizon for our metric is 
given by 



-2 



cdt d(l+7) d(i+7>- 

<> - ct (25) 



tdd+^D d{l + i) 
so the equation may be recast as 

E = Eq(x, t) e l{k - r - k - Ld) = E (x, t) e ik ^ kx - L) (26) 
exactly similar to the Newtonian result where the horizon is simply L = ct. 

4. Electrostatic Oscillation 

In this section we shall very briefly consider an electromagnetic wave in a 2-component 
plasma. For simplicity we assume a small amplitude wave such that, v x B = be- 
cause it is of second order in perturbed quantities. This, in turn, allows us to neglect 
the motion of ions. Skipping intermediate mathematical steps for space we get via 
equations (14) and (16) for our metric (5) 



^ d(1+7) E = —^E (27) 



t ( not to scale ) 



Figure 2: oj d ~ t graph. As the number of dimension increases red shift decreases 



The last equation is very similar to the flat space case except that here, uid is not a 
constant but shares the background expansion. When one uses the last equation in 
Maxwell equation (10) through J = n Q qv we arrive at Coulomb's law 



VxB= -— eE (28) 
c 



where 



e(u d ) = (29) 



CO 



is the dielectric constant of the plasma medium with the suffix 'T' signifying trans- 
verse mode. The u 2 T is related to the well known plasma frequency [23], 



u 2 



b d n q 2 2 

(d 



(d-l)\\ 

/d+l) (d-i) 
2 y 2 7T 2 

(d-2)!! 



(d odd) 



Skipping mathematical details it can also be shown that the transverse and the 
longitudinal modes give the dispersion relations as 

0,2 = u 2 + c 2 k 2 ^ 

u 2 L = u 2 p (31) 

These relations are pretty well known in the special relativistic case, excepting that 
here all the quantities depend on time as well as the total number of dimensions. 



Apparently the equation (29) has the same Newtonian form but both the fre- 
quencies depend on the scale factor, l A' which, again, is a function of both the 
number of dimensions and the equation of state chosen. To end the section let us 
investigate the time dependence of the dielectric constant in equation(29). Now, 
LVp T should share the time evolution of the background electron number density, Uq 

I 2 

( the inverse of volume of the universe) i.e., uq ~ A" ~ t Again, from (27) 
we get, vT ~ A n<yl ~ d \ The equation (10) further dictates that v ~ A' 1 , which gives 

r ~ A (2 - d) . So, ~ 2a ~ A" 2 ~ t' 1 ^. On the other hand the equation (24) 
implies that uj\ ~ A~ 2 , hence e{ujd) does not explicitly depend on time. This is a 
remarkable result in the sense that for a FRW type of metric the dielectric constant 
is a real constant irrespective of not only the total number of dimensions but also 
on the equation of state i.e., e(u^) continues to remain constant all through the 
evolution of the universe. 

5. Electromagnetic Oscillations in Cold Plasma 

In this section we investigate the situation where a plasma in thermodynamic equi- 
librium is slightly disturbed through the passage of an electromagnetic wave. We 
assume that an external ambient magnetic field is also present. We, however, as- 
sume the plasma medium to be cold so that the pressure can be neglected when 
considering the particle equation of motion. In stellar systems one often encounters 
situations where relaxation times are much larger than the age of the universe so 
that collisions ( hence pressure ) may be neglected. The effect of an electric field is 
not generally seriously considered because of the well known Debye shielding effect. 
The general problem of an electromagnetic wave propagating along an arbitrary 
direction with the external magnetic field is given by Appleton and Hartee in the 
Newtonian case when studying the propagation of radio waves in ionosphere. Hol- 
comb [11] studied in FRW metric a specialised situation of the A-H equation in 
the dust case. Considering the fact that a general solution with arbitrary 9 is very 
difficult to tackle in an expanding background with arbitrary number of dimension 
we shall restrict ourselves to the cases when the electromagnetic wave propagates 
parallel and perpendicular to the magnetic field. However the topic is of great im- 
portance in astrophysics and space science where electromagnetic wave propagation 
in magnetized plasma is very relevant. 



we assume that the external, uniform magnetic field and the wave vector k are 
both aligned along the i th direction (say z direction with i = 3) in the d- dimensional 
space, such that k = \k\e 2 and B = |5|e z . As is customary in the analogous 3- 
dimensional static space we also assume that all the perturbed quantities have the 
same time dependence given by equation (23) such that the linearized equation of 
motion (13) takes the form 



[ Here E is _Lr to k and considering that we are dealing with a (d + 1) dimensional 
space time it has components E\, E 2 , E 3 , E±, . . . , Ej . ] 



Case I (B || k) : 




(32) 



Replacing ^ via equation (23) by 



dE 
~dt 



iuJd + 



E 



(33) 



which, when plugged in equation (10) gives, after a long but fairly straight 
forward calculation gives for j = l 



iu d 



A 



1-^ 

w2 



E x - 



<4 
^d - W ? 



—I- 



.UJd 



J'=2 

2 iA? 



(34) 

For the case j = 3 (i.e., along the direction of the magnetic field ) it takes a 
simple form 



(V x B) 3 = -i 



(35) 



repeating the process for the remaining (d — 2) components we can write for the fi th 
component a tensorial relation as 



(V x B)^ = -i^ v E u 



(36) 



(//, v = 1, 2, 3, . . . , d). where e^ u is rank 2 skew symmetric tensor of order 'd'. A 
little inspection shows that 



en 



£22 — ^44 — ^55 



£12 


= e i4 


= £i5 = . 


£31 




= £34 = • 


e 33 


= 1 - 


u 2 

— =P3 



edd = 1 



u 2 d - u 2 c 



Pi (say) 



Cld — o 9 ~~ ^2 



e 3 d = 



(37) 

(38) 
(39) 
(40) 



so the (d x c?) permitivity tensor comes out to be 

/ pi ip 2 ip 2 
—ip 2 pi ip 2 
p 3 

-ij3 2 -«P2 pi 



\ -iPd 

Here cu„ is the plasma frequency given by 



b d n e 2 







Pi ) 



(41) 



(42) 



and the electron cyclotron frequency is given by 




Figure 3: uj c ~ t graph for both radiation and matter dominated era. As the number of 
dimension increases uj c decreases. Further the decay is sharper in dust case compared to 
radiation case. 



eB 



eB 



m P c 



(43) 



where B, the orthogonal magnitude of the ambient magnetic field is given by B 

2 

| [B) Z {B) Z \ 1 / 2 = Bt d ^+^ for our system (see figure 3). 

If the magnetic field is switched off (u> c = 0) the equation (34) reduces to 



(V x B) j = -^eE 3 



(44) 



exactly similar to the expression(28) of the section 4. Thus the introduction of the 
magnetic field generates varied modes transforming the dielectric constant scalar e in 
equation (35) to a second rank tensor e^. Although the equations (30) - (40) exactly 
resemble the analogous expressions in Newtonian theory the fact remains that that 
all the frequencies now depend on time rather than being constant. Further the 

2(d-l) 

cyclotron frequency u c decays as t d ( 1 +T) exactly similar to the orthogonal component 
of the magnetic field. 

If we take the curl of the equation (36) and replace V by ike z , after a long but 
fairly straight forward calculation we are led to the matrix form 



/ (1 



El' 

r.2 , 



,'P2 

L n2 



P2 



Pi/ 
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pi > 
73 



( Ex\ 

E 3 



/0\ 






J \E d ) \ ) 



(45) 



where 



2 U2 



n 



c 2 k 



UJ- 



(46) 



where n is the refractive index of the plasma medium. 
Three modes are possible. 

First the longitudinal mode characterized by E\ = E 2 = E 4 = . . . — E d — 0, 
(£"3 ^ and pz = 0). Since the displacement is along Z direction the magnetic field 
has no role to play and u d = uo p . 

As we are more interested in the dynamics of the electromagnetic waves rather 
than the plasma oscillation as such we take E% = in equation (45). Setting the 
determinant of the resulting (d — 1) x (d — 1) matrix in equation (45) to zero we get 



n = p l ± p 2 
The plus sign gives via equation (36) 



(47) 



E^-iE u = (p^u^3 



such that 



Ei = (eft - ie„) 



i(k L z-u) d t) 



(48) 
(49) 



corresponding to left circularly polarized wave. 

The wave number ki can be found from equations (37, 38, 46, 47) as 



Kr, = 



to: 



1/2 



U) d (0J d + u c 

On the other hand for the minus sign in (47) we get 



E R = (e M + e v ) e [ 



OkRZ-Ut)} 1 / 2 



representing RCP wave with 



k R = — 



1/2 



(50) 



(51) 



(52) 



It is clear that the two eigen modes have different phase and group velocities and 
unlike the former case the Right Circularly Polarized (RCP) wave has a resonance 
at uj = ujf = uj c where the phase velocities vanish. The expressions so far exactly 
resemble the ones found in the propagation of an electromagnetic wave with an 
ambient magnetic field in Newtonian mechanics. However, here all the quantities 
Ud, oj p etc. depend on time, a dependence modelled by the form of line-element, the 
number of dimensions and also the equation of state. 

It should be noted that the time dependence of Ud and uj c are different being 

2 _2(d-l) 

ujd ~ t d ( 1+ T) and u c ~ t . So there is no fixed resonant frequency as in 
Newtonian case but with time it changes. 
It also follows from equation (52) that for 



UJl = - 

2 



(53) 



the wave vector k vanishes and our analysis breaks down. So the wave propagates for 
< uj c and io-y < Ud, otherwise it becomes evanescent. Moreover, as the temporal 

dependence of uj c and Ud are different the magnitude of the allowed region changes. 

With dimensions uo c decays more sharply than Ud- Thus the propagation of the 

electromagnetic wave is more restricted in higher dimensions than the usual 4D. 
Returning to the Left Circularly Polarized (LCP) wave we see that the wave 

vector vanishes for 



Ud = ^2 = 



(54) 



and so the wave propagates for ujd > uj 2 

From what has been discussed above it is tempting to look for Faraday rotation 
(see ref. 4 for recent astrophysical data) analogous to the Newtonian case. Assuming 
that an electromagnetic wave traverses a distance z in a plasma medium with a 
magnetic field subject to the restriction on frequencies discussed above the Faraday 
rotation is given by 

0= kL_kR z (55) 
It should be noted that one should revert to the physical co-ordinate rather than 

2 

the co-moving one we are considering here. Accordingly Z p ^ = t d (!+T) , Z cm and 8 
finally comes out to be via equation (23) 

9 = ^^Z ph (56) 

with no dependence on time, so apparently the number of dimensions and the equa- 
tion of state have no impact on this classical result. It may be relevant to mention 
that measurements of the radio waves from the extra galactic sources suggest that 
various spiral galaxies are endowed with magnetic fields whose intensities are of the 
same order of magnitude as that of Milky way [1] i.e., of the order of microgauss 
corresponding to an energy density stored today in CMBR energy spectrum peaked 
around a frequency of 30 GHz. 

Case II (B±k) : 

Let us very briefly consider the case of a plasma with a uniform magnetic field 
B = B e z , through which an electromagnetic wave is propagating with propagation 
vector k = ke x , perpendicular to the magnetic field. Here two modes are possi- 
ble. As the mathematical exercise closely resembles the case I we shall totally skip 
intermediate steps to write the final form as 

1. First mode (called ordinary wave) with displacements in z direction ( i.e. || B) 
having the dispersion relation 

= 4 + Pc 2 (57) 

as the magnetic field has no influence for motion parallel to itself the equation (57) 
is exactly same as equation (30) for the electrostatic oscillation. 



2. Second mode (called extraordinary wave) with displacements in (d-l)-dimensional 
hypersurface (-Li?) having dispersion relation 

= 1 _ ^ <4 - <4 (58) 

Before ending the section a final remark may be in order. We know that pulsars 
are rotating neutron stars giving out pulses of radio waves periodically, which are 
affected by the interstellar medium during their propagation to reach us. If the 
interstellar medium has a component of magnetic field parallel to propagation di- 
rection then as shown earlier the plane of polarization will suffer Faraday rotation 
depending on frequency, having a spread in the rotation angle. This spread may 
have, in principle at least, some imprint on the nature of expanding universe. 



6. Discussion 

With the help of (3+1) formalism the Einstein-Maxwell and the electrodynamical 
equations are written for a (d+1) dimensional FRW-like spacetime in presence of 
plasma and linearised equations are solved for different phases of the universe. The 
analysis essentially generalises to HD the well known results of Holcomb and Tajima. 
The salient features of our analysis may be summarised as: 

1. For a propagating wave in HD in vacuum the photons redshift most in 4D 
and for a fixed d in radiation dominated model. 

2. Although the plasma is sharing the expansion of the background the dielectric 
constant remains a true constant. So the photons are in a sense self similar. This 
result was found earlier by Holcomb and Tajima. We here generalise this remarkable 
result to the case of extra dimensional spacetime and also for a fluid obeying a general 
equation of state. It may be tempting to suggest that the fact that the classical flat 
space result of the constancy of the dielectric constant is carried over to non static 
curved background and that too in higher dimensions may be due to the conformal 
flatness of the particular metric analysed here. So one should move with caution 
against any far fetched generalisation and in other complicated space time this result 
may not be true. 

3. In the presence of an external magnetic field many interesting oscillation 
modes manifest themselves. A simplified Appleton-Hartee type of solution gener- 
alised to higher dimensions is obtained in curved spacetime. Only a selected range 
of frequencies are available for propagation here. 

4. The well known phenomenon of Faraday rotation is obtained. 

To end a final remark may be in order. The present work suffers from two se- 
rious disqualifications. For sake of mathematical simplicity we work out everything 
assuming a linearised plasma theory. Conditions under which one may assume lin- 
earized plasma theory may well exist in Newtonnian theory, but we are not being 
able to clearly formulate those things for the case of a nonflat spacetime and that too 
when it is expanding. Secondly most observational evidences suggest that even if 
one starts with a higher dimensional phase the universe underwent the self compact- 
ification transition much earlier than the epoch when the big bang nucleosynthesis 



sets in. So although literature abounds with works (for example, higher dimensional 
black holes and its thermodynamics etc.) studying the standard electromagnetic as 
well as MHD laws in the framework of multidimensions becomes a sort of suspect . 
In that sense our analysis is more of a purely theoretical nature without much direct 
physical applications. In future work one should try to generalise these results in 
the realm of non linear plasma and also attempt to relate some of our findings to 
known astrophysical data. 
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